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Abstract 

A non-string-theoretic calculation is presented of the microcanonical entropy of relic 
integer-spin Hawking radiation — at fixed total energy E — from an evanescent, neutral, 
non-rotating four-dimensional black hole. The only conserved macroscopic quantity is the 
total energy E which, for a black hole that evaporates completely, is the total energy of the 
relic radiation. Through a boundary-value approach, in which data for massless, integer- 
spin perturbations are set on initial and final space-like hypersurfaces, the statistical- 
mechanics problem becomes, in effect, a one-dimensional problem, with the 'volume' of 
the system determined by the real part of the time separation at spatial infinity - the 
variable conjugate to the total energy. We count the number of field configurations on the 
final space-like hyper surface that have total energy E, assuming that initial perturbations 
are weak. We find that the density of states resembles the well-known Cardy formula. The 
Bekenstein-Hawking entropy is recovered if the real part of the asymptotic time separation 
is of the order of the semi-classical black-hole life-time. We thereby obtain a statistical 
interpretation of black-hole entropy. Corrections to the microcanonical entropy are com- 
puted, and we find agreement with other approaches in terms of a logarithmic correction to 
the black-hole area law, which is universal (independent of black-hole parameters) . This 
result depends crucially upon the discreteness of the energy levels. We dicuss the similar- 
ities of our approach with the transition from the black-hole to the fundamental-string 
regime in the final stages of black-hole evaporation. In addition, we find that the squared 
coupling, g 2 , which regulates the transition from a black hole to a highly-excited string 
state, and vice versa, can be related to the angle, 5, in the complex-time plane, through 
which we continue analytically the time separation at spatial infinity. Thus, in this sce- 
nario, the strong-coupling regime corresponds to a Euclidean black hole, while the physical 
limit of a Lorentzian space-time (the limit as S — > 0+) corresponds to the weak-coupling 
regime. This resembles the transition of a black hole to a highly-excited string-like state 
which subsequently decays into massless particles, thereby avoiding the naked singularity. 

1. Introduction 

Bekenstein [1] was the first to associate an intrinsic entropy with a black hole, such 
that black-hole entropy measures one's ignorance about the black hole's internal state. The 
logarithm of the number of distinct internal configurations of the hole, when one studies 
the microcanonical ensemble, in which the macroscopic external parameters (mass, charge 
and angular momentum) are fixed, then determines the microcanonical entropy. These 
internal configurations are related to different possible pre-collapse configurations, which 
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result in the formation of a stationary black hole with the same external parameters. In 
this interpretation, the black-hole entropy is the logarithm of the number of distinct ways 
in which the black hole may be formed by infalling matter. 

In relating the dynamical degrees of freedom of a black hole to its entropy, it was 
suggested that these degrees of freedom refer to states of all fields which are located inside 
the black hole [2]. In particular, modes of fields inside the black hole in the vicinity of the 
horizon give the dominant contribution to the entropy. Averaging over the states located 
outside the hole then enables the statistical-mechanical entropy of the black hole to be 
determined. 

The main aim of our approach is to give the black-hole entropy a statistical-mechanicalj 
interpretation in terms of a counting of states. We proceed by considering a (quantum 
or classical) boundary-value problem, in which data for massless, integer-spin, wave-like 
perturbations are set on initial and final asymptotically-flat space-like hypersurfaces Ej 
and E p , separated by a large proper time T at infinity [3] . (Fermions will be considered 
in another paper.) In [3], this led to calculations of the quantum amplitude (not just the 
probability) to go from (say) an initial spherically-symmetric pre-collapse configuration 
to a prescribed final distribution of radiative fields after the black hole has evaporated 
completely. In order to make the classical boundary- value problem well-posed, one must 
rotate the asymptotic proper time T into the lower complex plane; the real-T hyperbolic 
boundary-value problem for wave-like field perturbations is not well-posed [4]. 

Below, it will become clear that the real part of the asymptotic proper time T between 
the initial and final hypersurfaces plays a major part in the computation of the entropy. 
Apart from the total ADM energy - - the only macroscopic quantity appearing which 
encodes the three-dimensional geometry of the initial and final hypersurfaces — the only 
other macroscopic quantity is the one-dimensional 'volume' of the system, here effectively 
given by T . The 'volume', being conjugate to the fixed total energy at infinity, assumes 
a large range of values. However, we shall not take the thermodynamic or infinite-volume 
limit. We show that, when the 'volume' T is of the order of the evaporation time of the 
black hole, and when the given (fixed) total energy, measured after complete evaporation 
of the black hole, is E , then the logarithm of the number of final radiation microstates is 
given by the usual expression for the initial entropy of the black hole, which semi-classically 
has the form A-nks M 2 (m p )~ 2 . This gives a statistical-mechanical interpretation for the 
black-hole entropy. 

The final state on the final space-like hypersurface E^ is labelled by a set of variables 
{E, bj} . Here, the {bj} denote the massless, integer-spin 'Fourier' amplitudes [see Eq.(1.6) 
below], j labels the quantum numbers associated with the relic radiation, and E is the 
fixed total energy of the relic radiation. For simplicity, consider the example of small 
perturbations e<j)^ , about a spherically-symmetric background <&(£, r) , of a massless real 
scalar field = <fr + e (p^ , where £ is a small expansion parameter. One expands the 
perturbation function (p^ in the form [3] 






1=0 m=-£ 



(1.1) 



2 



The background metric is taken in the form 



ds 2 = - e h(t ^ dt 2 + e a ^ r) dr 2 + r 2 (d8 2 + sin 2 6> d(j) 2 ) . (1.2) 
The perturbed scalar wave equation at late times leads to the (£ , m) mode equation 

(e^ h - a ^ 2 d r ) R £m - (d t ) 2 R lm - - (d t {a-b])(d t R im ) - V,{t,r)R im = 0, (1.3) 
where 

v tM = ^( <(< + 1) + ^) . (,4) 
Here, m(t , r) is defined by 

exp(-a(t,r)) = 1 - 2m & r) . (1.5) 

Near the final hypersurface E p , one can decompose the radial wave functions 
Ri m (t,r) harmonically with respect to time, and the resulting functions Rktir) can be 
normalised in a convenient way [3]. With this normalisation, when one works near the 
final hypersurface E p , one can write 

1 f°° 

= - E / dk h ^rn Rki(r) Y tm (Q) . (1.6) 

This defines the bj coefficients for spin s = ; thus, j denotes the collective indices 
klm . The condition that E is fixed corresponds physically to a non-equilibrium situ- 
ation, where there is an exchange of energy between the black hole and infinity, which 
results (after a long time-interval) in the disappearance of the black hole (and possible 
later re-appearance, etc.). As there are no other conserved quantities apart from the total 
energy, an isolated system in equilibrium, having total energy E , would sample all its 
eigenstates at that energy with equal probability. The amplitudes {bj} , however, corre- 
spond to different final field/particle configurations (microstates) but identical total energy 
E (macrostate). The coarse-graining arises from the loss of phase information about the 
final configuration. The microcanonical entropy, therefore, is the logarithm of the number 
of different microstates that correspond to the same macrostate. 

One would also like to understand in more detail what happens in the final stages 
of black-hole evaporation. For example, it is expected that, within a semi-classical and 
purely Lorentzian picture, a naked singularity would form just as the black hole finally 
disappears. Were the black hole to disappear completely, then there would be a breakdown 
in the unitary evolution of quantum mechanics [5]. That is, pure quantum states would 
evolve into mixed states when the information about black-hole states located inside the 
horizon is lost. Here, though, we follow the approach of [3], in which the proper time 
T at infinity is displaced into the lower complex half-plane, and so any singularity in 
the classical solution is avoided. This leads us to consider the microcanonical ensemble; 
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pure quantum states evolve into pure quantum states, and there is no information loss. 
In another paper, we give in more detail the resulting unitary description of quantum 
black-hole evolution [6]. 

In Sec. 2, we compute the total energy of the relic Hawking flux crossing the final hyper- 
surface, in the limit of a large total energy or mass. In Sec. 3, we calculate the microcanon- 
ical density of states as a function of total energy and 'volume', that is, of the asymptotic- 
time interval T. In Sec. 4, we summarise the theory of the black hole/fundamental string 
transition in the final stages of evaporation, in order to exhibit the similarities with our 
approach in a recent paper [7], a cosmological analogue of the black hole/string 

transition was described. The idea is that massive elementary particles behave like black 
holes at large coupling. In string theory, there are very massive states, and large numbers 
of states of a given mass. In this approach, the large degeneracy is regarded as the origin 
of the Bekenstein-Hawking entropy. There follows a discussion of the thermodynamics 
in Sec. 5. Sec. 5.1 is devoted to deriving the semi-classical Bekenstein-Hawking (B-H) en- 
tropy formula through constraints on the one-dimensional volume. In Sec. 5. 2, we compute 
the corrections, beyond the semi-classical order, to our micro canonical entropy formula, 
finding agreement with other approaches as to the universal negative leading-order cor- 
rection. Sec. 6 puts some of the ideas of Sees. 3 and 5 on a firmer basis, in terms of the 
black hole/string transition theory of Sec. 4. Here, we use the probabilistic interpretation 
of the theory, derived from the quantum amplitudes associated with the boundary data, 
and make contact with the Euclidean approach to black-hole evaporation. 

2. Energy Constraint 

The total energy, E, of the radiation fields on the final hypersurface Si? , is given by 



where 7i = e~ b T tt , with T ^ being the lowest-order perturbed energy-momentum 
tensor, which is quadratic in the field perturbations. There is no background matter 
contribution, in view of the final boundary condition that only relic Hawking radiation is 
present at late times. The gravitational background geometry is described by a metric 7^ , 
which is spherically symmetric, but time-dependent, being of the form given in Eq.(1.2). 
The background approximation is described more fully in [3]. At late times, the background 
metric can be accurately described by a Vaidya metric [8] . To simplify the exposition, since, 
in principle, fields of different spins s may be present, we shall consider here only a massless 
spin-0 perturbation, writing the perturbative scalar field as A^(x); the further inclusion 
of spin-1 and spin-2 perturbations is straightforward [3]. On integrating by parts, using 
the linearised massless scalar-field equation and employing the boundary condition of 
regularity at the origin, {r = 0}, together with use of the asymptotic behaviour 
r 2 A^ (<9 r A (1) ) ~ r" 1 (or faster) at large r , we find, for our perturbative scalar field: 



E 




(2.1) 



E 




(2.2) 
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The second term represents the non-linear coupling between a changing massless scalar 
field and a changing background geometry. A feature of this term is its resemblance to 
the action integral for A^(x), save for the |(d — b) term in the integrand, which arises 
from the time-variation of the background metric. This interaction term is analogous to 
the phase and area/number interaction discussed in [9], which is there regarded as leading 
to the Hawking effect. The phase term is represented by the action, while the |(d — b) 
term is effectively the black-hole energy width, which describes the rate of change of the 
mass m(t, r) inside a radius r at time t in the Vaidya-like region, where e~ a ~ e b ~ 
1 — 2Gm(t,r) jr ; this term is identified with the area. Since we are here only concerned 
with physics as measured on the final hypersurface T,p at late times T, the width term will 
be slowly varying and of order 0(m~ s ). If one removes it from the integrand and estimates 
the action term as being 0(m 2 ) (on dimensional grounds), one finds that the interaction 
term is 0(m~ l ). In [9], the interaction term was indeed shown to be 0(m _1 ), thus giving 
a small contribution for large m. Writing u for a typical mode frequency of the perturbed 
scalar field [10], one finds that, in the adiabatic approximation u ^> \ |d — b\ , the energy 
E decomposes approximately into a sum over modes. 

Including now also spin-1 and spin-2 particles, with their polarisation states, the total 
energy of the final radiation field, in the adiabatic approximation, is 



E = M<? = h £ £ £ £ r~duu \b s ^ mP \ 2 , (2.3) 

S =0,l,2 £=s P m=-£ { >' J0 

where M is the initial ADM mass. 

The amplitudes {b SLU e m p} relate to the relic radiation reaching X + , and can be ex- 
pressed in terms of final amplitudes associated with the space-like hypersurface T, F [10]. 
As described in [3], our boundary-value approach leads to a discrete set of final radiation 
frequencies, given in the adiabatic approximation by {u n = 2iinT~ 1 ; n = 0, 1,2, . . .}, 
where T = 2 |T| and where T is the (complex) asymptotic time-interval at spatial infinity, 
between the initial and final boundary hypersurfaces E/ and T,p ■ This gives 



oc 



E = n 1^ 1^ 1^ 1^ uTTsV 2s ~¥ Un \ bsnirnP \ 

s=0,l,2 e= s P m=-e V 7 ' n=l ± 



|2 



2-kK 
f 



(2.4) 



J2 n N n , 



n=l 



where 



w » = f E E £ E lw I 2 ^ ( 2 - 6 > 



is the occupation number of the mode n, and we have introduced a frequency cut-off 
through n max . Therefore, the total excitation energy of the gas of particles is given in 
terms of the total oscillation quantum number as 

E = E N = ™ , (2.6) 
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where 

max 

n=l 

Thus, several states have the same total oscillation number. This equation describes 
combinatorially the splitting of N into a sum of the numbers n = 1 , . . ., n max , with 
each number n appearing N n times in the sum. Incidentally, Eq.(2.6) also gives the mean 
one-dimensional density of states 

n(F) f{E) - V 

P{E) ~ — ~ 2^hc ' (2 - 8) 

where T is the period of the phase-space trajectory with fixed energy E and V = cT, 
defined below in Eq.(3.1). 

Naturally, the quantity n max is given by noting that the radiation energy should not 
exceed the total energy E. That is, 

^n max = E , (2.9) 
whence n max = N. Thus, Eq.(2.7) implies that the ensemble average value of N n is unity: 

N 

\ 



1 N 

- nN n = 1 • (2-10) 

n=l 

By ergodicity, the ensemble average is equal to the time average for large N. This 
result follows from Eqs.(2.6,10), on substituting iV for T, and then taking the continuum 
limit in the summation. 

If the relic Hawking-radiation state is independent of the azimuthal quantum number 
m, then one can write 

J s=0,l,2 P i=s 

where 

c se = (2t+l) , (2.12) 

allowing for the degeneracy in the angular momentum modes. With a view to simplifying 
the calculation, a further assumption would be to approximate the black hole as a point 
radiator, which emits only s-waves. S-wave emission is only sensitive to the radial coordi- 
nate, and higher angular modes are less significant. In general, a black body of radius R 
emits modes with angular momenta I < £ max = u> R, where u> is the radiation frequency. 
This would introduce a cut-off in the t summation. 

3. Density of States 
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As is well known, the canonical-ensemble approach to black-hole evaporation breaks 
down in Schwarzschild space-time, as the Hawking temperature, Tjj , is inversely related 
to the black-hole mass, M, whence the specific heat C = dM/dTn is negative. Another 
way of seeing this is to note that the density of black-hole states increases roughly at the 
rate e M , for large M, resulting in a canonical partition function which diverges at large 
M. Therefore, one should instead employ the micro canonical ensemble, which involves a 
fixed total energy E = M c 2 ; in this case, the total energy is not permitted to fluctuate. 
The micro canonical entropy does, however, involve corrections to the Bekenstein-Hawking 
area law, due to quantum space-time fluctuations at fixed horizon area (equivalent to fixed 
energy), as seen, for example, in the loop quantum gravity approach [11]. In Sec. 5. 2, we 
discuss the nature of these entropy corrections. 

In this Section, we calculate the microcanonical total density of states Q(E, V). For 
an isolated system, this counts the number of microscopic configurations in a volume 

V = cf , (3.1) 

possessing energies distributed between zero and E. The energy constraint is the only 
conservation law that we impose. In the quantum theory, when one is counting, different 
degrees of freedom correspond to distinct states. 

Another way of interpreting the microcanonical density of states is as 

oo 

Q(E,V) = Q n( E ' V ) > ( 3 - 2 ) 

N=l 

where Q.^(E, V), as below, gives the number density of states {N^} which contain N 
particles. These N particles may be labelled by frequencies, ordered such that 
< <^2 < ••• < • In Eq.(3.2), £l(E,V) gives the phase-space volume, up to 
the energy E , of the system with volume V , where N denotes the number of degrees 
of freedom. We are effectively considering a set of oscillators labelled by {bj}, whose 
phases are random; that is, their location in phase space is unknown or their phase-space 
trajectories are uncorrelated. If the phases were not random, the motion of the system 
would not be ergodic, and entropy could not be defined. 

In order to determine £l(E, V), we first define 

X(E,V) = sum {Nn} 0(E - hJ2 ^iV n ) , (3.3) 

n 

where 9 is the Heaviside step function and one sums over all states of the fields, labelled 
by {N^}. The density of states is defined to be 



n(E,V) = §se 



5E 5 { E ~ ft^niVn) 



(3.4) 



{N n } 
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namely, the total number of configurations in the energy range [E , E + SE] . The factor 
8E is some spread in energy or tolerance - independent of E, so as to make Q(E, V) di- 
mensionless that depends on the precision with which the system has been prepared. 
One apparently natural condition is that E ^> SE, and that, for any energy fluctua- 
tion, AE 3> SE. For discrete systems, however, SE has a minimum value, although for 
continuous systems it can be arbitrarily small — see Sec. 5. 2. 

The number density Q^(E, V) may similarly be represented in the form: 

Sl#(E,V) = SE S ( E ~ h J2 ^niVn) • (3.5) 

{N N ,} V n=l / 



This expression (3.5) may also, in the continuum limit, be re-written as: 

se r E * f \ 

n#(E,V) = -T- / H dE n p(E n ) slE - J2 E n) , 

^ ■ ^° n=l ^ n=l ' 



(3.6) 



where E n = hu) n N n , and where p(E n ) (the single-particle degeneracy) is given by 
Eq.(2.8). The factor Yln=i dE n p(E n ) counts the number of iV-particle states such that 
the n-th particle has energy E n . The l/(iV!) term accounts for identical particles in 
classical statistics, while the delta-function constraint counts the number of iV-particle 
states with energy E ; for further details, see [12]. 

Thus, we consider the set of states {A^ n } as comprising subsets, labelled by N, such 
that Qfi(E,V) is the number density of states possessing N particles. One can then 
evaluate Q^(E, V), giving 



AT-l 

SE\ 1 (T E 



E 1 
n#(E,V) = SE —, II — 



E J N\{N-1)\ V h 
Now, define a (small) dimensionless parameter g by: 



N 



(3.7) 



9 = — = N ■ (3 ' 8) 

where the second equality comes from using Eqs.(2.6, 3.1). We then obtain the energy 
spectrum [13]: 

E = gE p VN , (3.9) 

i 

where E p = (hc 5 /G) 2 is the Planck energy. Notice that, for E 3> E p , the density 
£lfi(E,V) increases as g 2 decreases, corresponding to the notion that the phase-space 
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volume and therefore also the entropy increases as the observation time T 

exceeds the dynamical time-scale of the collapsing matter, as the black hole forms and 
loses its hair. In a related calculation in a cosmological context in [14], the aim was to 
find a gravitational entropy function which increases monotonically as the system becomes 
more inhomogeneous. 

The microcanonical entropy is then given by: 

S(E,V) = k B logf n^(E,V)) . (3.10) 

One can approximate the right-hand side of Eq.(3.10) as 

S(E,V) ~ k B log(fi#(E,V0) , (3.11) 

where iV is the value of iV which maximises log(o,^(E, V)^; that is, the most probable 
value of N : 

8fl - 

-f = N(E,V) = . (3.12) 
oN n 

Using Stirling's formula for large N, namely, N\ ~ y/2n e~ N N N+ ^ , and assuming that 
5E/E is independent of A^ (see Sec. 5. 2), we obtain 



/ EV 

N * yj— » 1 . (3.13) 

Note that N = 1 corresponds to an elementary particle, as T is then the black-hole 
Compton time. Eq.(3.13) seems to be related to the proposal in [15] for the maximum 
number of space-time bits in a spatial region with total volume V and energy E. Therefore, 
the degeneracy is given by 



n(t\V) - cxp(2 A /^) , (3.11) 



and 



EV 

S(E,V) ~ 2k B \l— = 2k B N = 2h B \fN . (3.15) 



This is one of the main results of this paper. Note that the entropy is independent of the 
gravitational constant G . 

In arriving at Eq.(3.15), the 5E/E term was neglected in comparison with the large 
N term. This will be justified in Sec. 5. 2. Thus, we find the well-known result that, for 
ideal bosons in a one-dimensional harmonic-oscillator potential, the asymptotic density 
of states, as A^ — > oo , is equal to the number of ways of partitioning an integer A^ into a 
sum of other integers [16]. 
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Through a little manipulation, the entropy may be expressed in the form 



S = 



2k B 
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(3.16) 



Thus, the entropy is linear in the total energy E . The (7-dependence illustrates that, as 
g — > 0, the entropy increases, as one expects in the classical limit. 

4. Black— Hole/String Correspondence 

One could compare the results in Eqs.(3.9,15) with those in [17], which assert that 
the entropy of a Schwarzschild black hole in string theory is proportional to , where 
N s is the level number of a long excited string which collapses into a Schwarzschild black 
hole at a critical string coupling. 

Susskind [18] originally developed the idea that there is a one-to-one correspondence 
between quantum Schwarzschild black holes and fundamental string states, building upon 
the ideas first developed in [19]. One immediate similarity is that the degeneracy of states 
for a given mass increases very rapidly with mass, both for black holes and for fundamental 
strings. Yet, the leading string entropy term depends linearly upon the mass, whereas the 
Bekenstein-Hawking entropy is quadratic. The Correspondence Principle [17], however, 
aims to match these differing behaviours for a particular string coupling g c . A related 
(semi-classical) quantum field theory /quantum string black hole duality was investigated 



The one-dimensional string is characterised by its tension ~ (£ s )~ 2 , where £ s is the 
string scale, and by a (dimensionless) interaction strength g . Varying the coupling g has 
an interpretation in terms of a varying dilaton background, where the dilaton is given by 
a scalar field (/>, through the relation g 2 = e^. In practice, we treat 4> as a constant and g 
as a parameter. Consider a highly-excited string state and a weak string coupling j < 1, 
where, in four dimensions, £ s and the Planck scale £ p = ^/Gh/c 3 are related through 



Thus, we initially consider a string scale far above the Planck scale. By increasing 
the coupling adiabatically, (in effect, then, by increasing the strength of the gravitational 
interaction), one reaches a point such that the Schwarzschild radius of the string, given 
in terms of the string mass scale M s by 2GM s /c 2 , becomes larger than £ s ; at this point, 
the string makes a (smooth) transition into a black hole. Conversely, however, had one 
started with a black hole, and then decreased the coupling, one would find that the size 
of the black hole would eventually fall below the string scale. The correspondence point 
would be that at which the black hole makes the transition to a string. 

The (closed) string mass spectrum has the form [17] 



in [20]. 



M 




(4.2) 
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where m p = h/ (£ p c) is the Planck mass, and N s 3> 1 denotes the level number of a 
long (highly-excited) string. Therefore, for fixed N s , the string mass increases as N s 
increases. The string spectrum involves an infinite tower of massive states, as well as 
gravitons and massless spin-one particles (for example). For states in free string theory, 
with a given mass, the most numerous states are the single-string states, since the string 
entropy depends linearly upon energy. The resemblance of Eq.(4.2) to Eq.(3.9) is evident. 

The Bekenstein-Hawking entropy of a four-dimensional Schwarzschild black hole, 
Sb-h , is given in terms of t v as 

Sb __ h = 4 ^; 2 fa> v , ,4.3) 

where M is the black-hole mass. Thus, at the correspondence point where M = M s (we 
assume that the mass does not change over the transition), Eq.(4.3) becomes 

4irk B (£ P ) 2 N S 

S B — h ~ = 4nk B g 2 N s . (4.4) 

(A) 

Therefore, at the critical coupling [17], 

(g c ) 2 ~ (N a )-* , (4.5) 

the Bekenstein-Hawking entropy becomes the string entropy 

^string ~ ^k B y/Nl , (4.6) 

which resembles Eq.(3.15). One sees that the degrees of freedom of an excited free string 
can match the entropy of a four-dimensional Schwarzschild black hole. It can also be 
demonstrated that the sub-leading entropy terms, which depend upon the mass, also 
agree [21]. In this case, the degrees of freedom of a perturbatively interacting string can 
be identified with Schwarzschild black-hole states. 

Instead of fixing g , and letting the black-hole mass M (and therefore N s ) 
vary as the black hole evaporates, one can, equivalently, keep the entropy fixed, 

such that Spring = Sb-h (equivalent to keeping the state fixed), while varying 

_i 

g adiabatically [22]. Note that, for g 2 < (N s ) 2 , one has Spring > Sb-h , and 

_ i. 

conversely, for the case g 2 > (iV s ) 2 , one has Spring < Sb-h ■ Thus, if there is 
some critical mass, M c , such that Spring > Sb-h , then a black hole can make a 
transition to a higher-entropy string state. This mass M c satisfies GM c /c 2 ~ £ s , 
that is, 

M c = M s g~ 2 . (4.7) 

This mass scale is much larger than the Planck mass m p , in the case that j € 1. Further, 
at this mass, the Hawking temperature equals the Hagedorn temperature 

T s = (hc)/(k B be s ) , (4.8) 
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where b is a constant which depends on the particular model and on the dimension taken 

for space-time [18,19]. Dimopoulos and Emparan [23] found that a black hole makes a 

transition to a state of highly-excited and jagged strings (a string ball) at the mass M c . 

The production cross sections of black holes and string balls can be matched at M c . The 

string ball then loses mass through evaporation, and at a mass just below M c , it 'puffs up' 

i 

to a larger random-walk size ~ i s (M/M s ) 2 (see below [23]). The string ball continues 
to evaporate with a quasi-thermal spectrum, at the Hagedorn temperature, during which, 
its size reduces to ~ £ s . The black hole can be approximated as a point radiator, as its 
wavelength is of order (he) / {k-B T s ), which is larger than the size of the black hole. Since 
the specific heat of massive string modes is negative, these cannot remain in equilibrium 
with an infinite heat bath. Therefore, they evaporate into massless modes, so avoiding the 
possibility of a naked singularity [19]. 

In [19], the assumption was made that there is an intermediate heavy-string state, 
which decays at the Hagedorn temperature into massless strings; that is, into elementary 
particles such as photons. Zalewski [24], however, suggested that the black hole transforms 
directly into massless string states, using statistical arguments and the property that the 
evanescent black hole should have zero angular momentum. 

In fact, one can also match the total width for string decay (into the dilaton, graviton 
and massless fields) with that for quantum black holes [25]. The total string width has the 
form [26]: 

T s = Xg 2 M , (4.9) 

where A is a dimensionless numerical factor. A formula analogous to Eq.(4.9) encapsulates 
the decay widths of all the Standard-Model particles, as well as of topological and non- 
topological solitons, and of cosmic strings [25]. The black-hole decay rate has the form 



BH 



1 dM 



M dt 

For the low-energy canonical ensemble, we use the Stefan-Boltzmann law 

dM 4tt a (R s ) 2 {k B T BH ) 4 

dt ~ c 4 h 2 



(4.10) 



(4.11) 



where a is a dimensionless constant, the Schwarzschild radius is denoted by 
R s = 2GM/c 2 , and T BH = (he 3 ) / (8irGMk B ) denotes the Hawking temperature. Thus, 

Tbh = {k B T BH f ■ (4.12) 

As T B h — ► T 8 and R s — > £ s ; that is, as the black hole's temperature increases towards 
the string temperature during evaporation, then 

Tbh - ^9 2 M S , (4.13) 
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which resembles Eq.(4.9) in the case that M ~ M s . 

A random-walk interpretation of strings at zero coupling also exists [22,27]. From 
Eq.(4.6), with n = y/N^ , the string entropy can be represented as an integer n, up to 
a factor, giving the number of random-walk steps or number of string bits in a polymer 
string representation [28] . Each step has length t s and mass M s . At zero coupling, the free 
string effectively has a random-walk size y/n £ s = (N s )* £ s 3> £ s , and the total length of 
the string is L = n£ s . 

5. Thermodynamics 

Evidently, from Eq.(3.15), S{E, V) is an extensive function of E and V, both of which 
determine the global properties of the system at infinity. The corresponding thermody- 
namical quantities, such as temperature 9, pressure P, specific heat Cy and average total 
particle number N, can be computed using the standard formulae 



where (3 = giving 



d 2 S \ 

dE 2 



v 



N = 



(3PV 
k B 

(5.1) 



6 = = T^rv/iV, P = ^ = k BP , C V = S(E,V)>0, N = N , 

Kb V V ksV V 

(5 ' 2) 

where p is the energy density. The second equation in Eq.(5.1) describes an incompressible 
fluid, for which sound propagates at the speed of light the Zel'dovich equation of 

state. We consider both E and V to be large, such that p is finite and E ^> hc/V . 
The thermodynamic limit, in which V is regarded as being arbitrarily large, should not 
be taken, as our approach is based on considering only finite time-intervals at infinity. 
Hence, the ideal gas approximation holds, in that interactions between particles are being 
neglected: 

PV = Nk B 9 . (5.3) 

The number N = E/(k B 0) is also the average number of particles with energy k B 6 in 
which the black hole of total energy E decays. The property that Cy > implies that a 
stable thermal equilibrium can be reached by the radiation at any temperature. 

One can express the entropy in terms of the temperature 6 and the volume or length 
V, as (see [29]): 



{ks? 

TV He 



9 V . (5.4) 



Thus, the entropy is proportional to the length V, and one may interpret the entropy as 
comprising bits of length V. One obtains the familiar expression for the entropy if one 
makes the replacement 

V - \^ 2 fL , (5.5) 
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where L is a length scale and / denotes the number of field species. Thus, 




where p = (EL) / (he) is dimensionless. This is the Cardy formula for central charge / 
[16]. Eq.(5.6) would have been obtained if, instead of summing over {N n } configurations, 
we had summed over {N^} configurations such that N n = X^=i ■ 

When the particle energy hui n is much larger than the corresponding temperature 9, 
then, for E ^> Hu> n , one has 

N > n > y/N , (5.7) 



and the microcanonical ensemble should be used. Further, in the very-high-energy tail 
of the spectrum, where quanta have energies comparable with the total energy E of the 
system, one has 

n ~ N . (5.8) 
On the other hand, for low energies which obey the condition 

n < < N , (5.9) 

the canonical ensemble description of the gas is valid [30] . 

Another interesting relationship, following from Eq.(5.2), is 

= hc(3y/N = g- x t v \TN , (5.10) 

which relates the wave period to the inverse temperature. Thus, in this representation, one 
interprets V effectively the time-separation at infinity as a root-mean-square 
distance, in which case it represents the distance covered by a Brownian random walk with 
step-size g~ l £ p , taking iV to be the total number of steps [31]. The analogue of the string 
length, i s , in our theory, is provided by i s , where 

L = ^-N~^ = hc/3 , (5.11) 

which is just the random-walk step size. In addition, with i s = H/ (M s c) , the mass scale 
M s , analogous to the string mass scale M s , is 

M 

M s = —= = gm v . (5.12) 



From Eq.(5.2), one has 

9 = TT ' (5 - 13) 

k B t s 
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which is the analogue of the Hagedorn temperature given in Eq.(4.8). 

In the light of Eq.(5.11), for the Euclidean quantum gravity approach to black-hole 
evaporation, the imaginary-time period is given by the inverse temperature of fields on 
the background space-time [32]. 

5.1 Bekenstein-Hawking Entropy 

Given a four-dimensional, non-rotating, neutral black hole with horizon area A, the 
Bekenstein-Hawking entropy has the well-known form 



S B __ H = ^4 = tetaC-V (5.14) 
4{i p f 



This is a semi-classical result, valid for A ^> (i p ) 2 ■ Let us now equate Eq.(3.15) with 
Eq.(5.14), which is valid for M > m p : 

S(E,V) = S B — H . (5.15) 

Matching the four-dimensional entropy with a two-dimensional entropy is just one way 
of comparing the thermodynamics of any two physical systems. Strong matching might 
involve equating partition functions which, then, leads to the matching of all thermodynam- 
ical quantities [33]. The weak matching condition, Eq.(5.15), which does not necessarily 
imply the matching of other thermodynamical quantities, implies that 

% G 2 M 3 

T(M) = (2tt) 3 ^-^ . (5.16) 

In the semi-classical approximation, the evaporation time-scale, to , of a non-rotating, 
neutral black hole of initial mass M, is [34] 

to = (3 7 ) -j^r > ( 5 - 17 ) 

where 7 is a dimensionless coefficient, which depends on the number of particle species 
which can be emitted at a given temperature, and on grey-body factors. Thus, T is of the 
order of to • 

Another interesting way to arrive at a time-separation at infinity of the order of to > 
comes through setting, in Eq.(5.10): 

^ - ^7W « 1 ■ (5 ' 18) 



in which case we arrive back at Eq.(5.16). In other words, the B-H entropy equals the 
1-dimensional entropy Eq.(3.15), when Eq.(5.18) is satisfied. Therefore, we seem to have 
enough configurations on to reproduce the entropy of black holes. However, if one 
considered the case in which 2tt g 2 VW ^> 1 , then one would still be in the black-hole 
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phase, since T < t . When 2tt g 2 y/N < 1 , one is in a post-evaporation phase (perhaps 
string-like). Eq.(5.18) is essentially the critical coupling condition of Eq.(4.5). 

The microcanonical entropy in Eq.(3.15) exceeds the Bekenstein-Hawking value, if T 
exceeds the right-hand side of Eq.(5.16): 

y > (^) 3 w- 2 • 

However, for T given by Eq.(5.16), one has, from Eq.(5.2), 9 = AT H , where T H is the 
Hawking temperature. 

We conclude that, if our final space-like hypersurface, Hp, were located so that the 
time separation at spatial infinity between the initial and final hypersurfaces were given 
by \T in Eq.(5.16), then the thermal entropy of the remnant radiation would be of the 
order of the initial black-hole entropy. In that case, therefore, the thermodynamics of the 
four-dimensional black hole would be encoded in that of a one-dimensional gas of non- 
interacting bosonic particles. There have been similar approaches that modelled black 
holes in any number of dimensions, using a one-dimensional gas of massless particles [35]. 

As an aside, we mention here briefly the connection between our approach and the 
idea of strings growing or spreading upon reaching the black-hole stretched horizon, which 
is located at a distance ~ l s from the event horizon [36]. One can re-write Eq.(3.8) in the 
following suggestive way: 

^ = cfT . (5.20) 

Setting tcross = 2r s /c , which is the time-scale for light crossing a Schwarzschild black hole 
of radius r s = 2GM/c 2 , then 

Across = 9 2 (^p) ■ (5.21) 

Due to the redshift effect, the string and all the information it carries appears, from the 
perspective of a distant observer, to grow as the string passes through the stretched horizon, 
while, to a freely-infalling observer, it appears to be a Planck-sized object as it crosses 
the stretched horizon. The string will spread across the entire black hole in a time [36] 

t a - ( l A 2 — - (5 22) 

where g is here defined as in Eq.(3.17). The time t spr ead is short, compared to the black-hole 
evaporation time to i provided that g is small. One can immediately see the similarities be- 
tween Eqs.(5.21) and (5.22), in the case that T is of the order of the black-hole evaporation 
time. Furthermore, both t C ross arid Spread seem to encapsulate the notion of information 
traversing the entire black hole. 

5.2 Corrections to the Microcanonical Entropy 
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Corrections to the B-H area law are important, and related to the evaporation of 
the black hole. These corrections to the microcanonical entropy, Smc , due to quantum 
fluctuations at fixed horizon area A , for A ^> (£ p ) , seem to have a universal form, that 
is, a form which is independent of the particular parameters of the black hole [37]: 

Smc = S B -h - g ^og(s B -H^ + constant + 0^(S B -h) , (5-23) 

where S B h is given in Eq.(5.14). Such corrections can also be interpreted as being due to 

thermal fluctuations about the black hole's equilibrium configuration. The — | log (5 b-h) 
term, in particular, holds also for all black holes whose microscopic degrees of freedom are 
described by an underlying conformal field theory, including BTZ and string-theoretic 
black holes [37,38]. This term also appears when one computes a corrected Cardy formula 
for the density of states [39] — see also [41]. 

We now compute the corrections to Eq.(3.15), of the type given in Eq.(5.23), which 
arise from quantum fields on the black-hole background (as studied above), and from the 
reduction in the area or mass of the black hole. The B-H law Eq.(3.15) is, of course, only 
valid for A ^> (£ p ) . One may either convert the N sum in Eq.(3.10) into an integral, and 
then perform a stationary-phase approximation, or one may compute the sum exactly. A 
straightforward calculation yields 



s(e,v) ~ 2kB y^ + \ kB log (V j^j + kBlog yw) + constant ' 

(5.24) 

where we have neglected higher-order terms. The + \ coefficient of the logarithm in 
Eq.(5.24) is in contrast to the results in [42]. It has been suggested that different quantum 
theories lead to different quantum corrections to the black-hole entropy [37] . 

One can, however, make further progress with the \og(8E/E) term in Eq.(5.24). Since 
the energy levels are discrete, SE has a minimum value, which one may take, naturally, to 
be the spacing between the particle energy levels: 

h (uj n+1 - 0J n ) = i . (5.25) 
Thus, for V > (hc)/(4E) and for 

he 

SE = hiv = — , (5.26) 

(that is, the condition that the energy spread depends solely on the volume), we find, up 
to unimportant constant terms, that 




(5.27) 
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which implies that 

S(E,V) ~ k B (2y/N- | log(2"/jv)^ + constant . (5.28) 

In terms of the leading — | log( ) term, this agrees with [37,39], provided that we identify 
Eq.(3.15) with the black-hole entropy. Clearly, the discreteness of the energy levels is a 
crucial element in deriving Eq.(5.28). 

One might compare Eq.(5.28) with the calculation of the logarithmic corrections to 
the Cardy formula for the density of states in a two-dimensional conformal field theory 
[39]. In [39], there is an arbitrary period, which is chosen so that the central charge is 
universal (independent of black-hole area); this again leads to the — | log(SB--ii) 
term - see also [40]. A similar conclusion holds in the formulation of the present paper, 
where the spread in energy is fixed in relation to the period at infinity between the initial 
and final space-like hypersurfaces, £/ and T, F , as in Eq.(5.26). 

6. Probability Distribution 

We now discuss further our previous results for quantum amplitudes in black-hole 
evaporation [3] . We set approximately spherically-symmetric boundary data on an initial 
space-like hypersurface, Ej , and weak, non-spherical boundary data on a final hypersur- 
face, E p , separated from E/ by a proper time T, as measured at spatial infinity. The 
initial data may be regarded as describing a nearly-spherical configuration of gravity and 
matter, shortly before collapse to a black hole. The weak final data may be regarded as 
describing a particular choice of very-late-time radiation. The quantum amplitude to go 
from the prescribed initial data on Ej to the prescribed final data on E^ is essentially 
given by expUSd), where is the action of a solution of the classical field equations 
which proceeds from the initial to the final data (if such exists). In particular, this is 
expected to hold in the case that the theory is invariant under local supersymmetry [3]. 
In the present case, this semi-classical expression for the quantum amplitude will only be 
used in the weak-field limit, as above [43]. In fact, the above boundary- value problem, of 
finding a classical solution to the wave-like (hyperbolic) equations of a field theory, is not 
well-posed [4] . The natural cure, following Feynman [43] , is to rotate the time-interval T, 
as measured at spatial infinity, slightly into the lower complex half-plane: 

T = \T\e~ iS , (6.1) 

where 5 (0 < 5 <C 7r/2) is a real phase. One finds that S c \ has both its real and 
imaginary parts non-zero. Hence, from the squared modulus of the quantum amplitude, 
one arrives at a probability distribution for the weak-field radiation, as measured on the 
final hypersurface E_p [44]. 

Consider a particular configuration of the / oscillators N„ (assumed to be indepen- 
dent), at each level n, associated with the boundary data on Ei?. As above, we write: 

/ 

N = N[K] , N n = , (6.2) 

i=i 
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The probability density has the form (cf. [44]): 




P[{K}} = Z f {g*)' exp --VV 7n ArM , (6.3) 



where 



7n = 2ng 2 n + 2e 2 p n , g 2 = - (6.4) 



and the 'partition function' is 



( ^ max y \ 

- E E ^ N n) ■ ( 6 - 6 ) 
n=l i=l / 

The parameter e is the small perturbation-expansion parameter given just before Eq.(l.l). 
The 27r g 2 n term in 7n arises from the contribution to the classical action from the timelike 
boundary E°° at spatial infinity, which joins Ej to E^ [32]. The e 2 term in 7n arises from 
the weak, non-spherical boundary data for integer-spin perturbations on Ej? [3]. 

The {N^} sum in Eq.(6.6) is over all sequences of independent occupation numbers 
, such that n > 1 , 1 < i < f . Despite the formal similarity of Eq.(6.6) to a thermo- 
dynamical partition function, one should view the parameter 2-7T g 2 as being conjugate to 
N, and not the energy. With Eqs.(6.1) and (6.5) in mind, this is analogous to a number- 
phase conjugacy in quantum mechanics, although 5 is not an oscillator phase. A standard 
calculation gives 

max £ 

Zf{9 2 ) = II (l-e" 7 ")" , (6.7) 

n=l 

and from Eqs. (6.2,3), we obtain the average 

N n = ~~~~ T ■ (6.8) 

e 7 ™ — 1 

A physical interpretation for g 2 (and hence for the phase 5) can be obtained by 
evaluating 

^ max 

N = E n *n (6-9) 

n=l 

and using Eq.(6.8). Assuming that one can convert the sum to an integral (which is the 
case when g 2 <C 1 ), we find, from Eq.(6.3), that 



2ng 2 ~ AT* . (6.10) 
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Alternatively, let us assume that A n > e 2 p n , where A = 2tc g 2 . Then, 



p[{K}} = z f (\) 



Define 



Then, 



p(N) 



Z f (X) 



s[n - N({K}) 



dN p(N) e 



— \N 



JO 



which is just the Laplace transform of p(N) . Therefore, 



(6.11) 



(6.12) 



(6.13) 



satisfies 



Let us make the Ansatz 



P(N) = Z f (X) 



dN P(N) 



XN 



p(N) 



p(N) = c N-* 13 exp^av^) 



(6.14) 



(6.15) 



(6.16) 



where cq , a and (3 > are all constants. Then, the probability density Eq.(6.14) is peaked 
around N = N , where 



N = 



1 

4A 



a ± v 7 " 2 - 8/3A 



(6.17) 



If we further assume that a 2 3> 8/3 A and that A^ 3> 1, we obtain 

a 



A = 2ng 



?,2 



2JN { 



(6.18) 



o 



This agrees with Eq.(6.10), provided that a = 27Ta///6 and N = N . The inequality 
a 2 3> 8/3 A seems reasonable for large /. Therefore, 



p(N) 



co 



N~*P exp ( 2tt 



6 



(6.19) 



The logarithm of Eq.(6.19) agrees with Eq.(5.6) to leading order. Eqs.(6.10,18) are analo- 
gous to the critical string coupling condition of Eq.(4.5). This suggests that the parameter 
g behaves like a string coupling parameter. 
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There are further indications that the rotation angle 5 is analogous to the (squared) 
string coupling. In the squeezed-states approach to black-hole evaporation [45], the 
squeeze parameter r,- is related to the frequency Uj and the black-hole mass M through 



tanhrj = exp ( - — 3 - ) . (6.20) 



Since < Nj >= sinh 2 (rj), this gives a thermal distribution 

< Ni > = -s—^ , (6.21) 

at the Hawking temperature Tjj = (ks Ph) 1 ■ This spectrum is only valid in the phase 
for which the black hole emits particles at low energy, (which, of course, corresponds to 
the majority of the lifetime of the black hole). For the final stages of evaporation, the 
energy spectrum does not have the canonical form, but a micro-canonical distribution 
[12]. In the approach of [44], the boundary-value problem of [3] for computing quantum 
amplitudes for wave-like massless perturbations at late times was related to the squeezed- 
states formalism. It was shown that the limit 5^0 for the rotation angle in Eq.(6.1) 
corresponds to a highly-squeezed final state for the relic Hawking radiation, with a squeeze 
parameter rj given by 

tanhrj = exp^— 2ujj \T\ sin^j . (6.22) 

Eqs.(6.20) and (6.22) agree when 

hN sin S 4nGM 



Mc 2 " c 3 



(6.23) 



Here, N is given in Eq.(2.6). For small 6, this gives Eq.(3.9), with g replaced with g. 

In addition, Damour and Ruffini's approach to black-hole evaporation [46] involved 
complexifying the Schwarzschild mass M, in order to smooth out large oscillations in the 
outgoing particle wave-function as one approaches the future event horizon. Whether 
one continues M into the lower or the upper half of the complex plane depends upon 
whether one wishes to describe positive- or negative-energy states, respectively. In this 
procedure, quantum states are smooth across the event horizon, and describe a negative- 
energy component that falls towards the curvature singularity, together with a positive- 
energy component that can escape to infinity (depending on its energy). If, following this 
approach, one sets 

M = \M\ exp(-id) , (6.24) 
then, for small 5, one has M = |M| — (i T)/2 , where we define 

T = 4g 2 \M\ . (6.25) 
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This expression bears a resemblance to the expression in Eq.(4.9) for the width or decay 
rate, taking g as in Eq.(6.5). 

This analysis suggests that the angle 6, by which the time-interval at spatial infinity, 
T, is rotated into the complex plane, plays a role analogous to that of the string coupling 
g. Black-hole evaporation can therefore be viewed in the following way: Initially, S is 
non-zero, and is given by its value in the Euclidean regime. In the limit 5^0+ (and, 
hence, for weak coupling g — > 0) , one moves towards the physical, Lorentzian, regime, and 
the black hole makes a transition to a 'string-like' state, where it subsequently evaporates 
into massless particles. 



7. Conclusion 

In this paper, we have expressed the black-hole entropy via a counting of massless 
field configurations which have a fixed total energy on a post-evaporation space-like hyper- 
surface. We have also explored some of the analogies with the string/black-hole correspon- 
dence, incorporating some of the notions of black-hole evaporation in a complexified-time 
representation. 
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